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background: Although the role of human papilloma virus (HPV) in cervical squamous cell carcinoma (CSCC) is well established, the 
role in head and neck SCC (HNSCC) is less clear. MicroRNAs (miRNAs) have a role in the cancer development, and HPV status may 
affect the miRNA expression pattern in HNSCC. To explore the influence of HPV in HNSCC, we made a comparative miRNA 
profile of HPV-positive (HPV+) and HPV-negative (HPV-) HNSCC against CSCC. 

methods: Fresh frozen and laser microdissected-paraffin-embedded samples obtained from patients with HPV+/HPV— HNSCC, 
CSCC and controls were used for microarray analysis. Differentially expressed miRNAs in the HPV+ and HPV— HNSCC samples 
were compared with the differentially expressed miRNAs in the CSCC samples. 

results: Human papilloma virus positive ( + ) HNSCC had a distinct miRNA profile compared with HPV— HNSCC. Significantly 
more similarity was seen between HPV+ HNSCC and CSCC than HPV— and CSCC. A set of HPV core miRNAs were identified. 
Of these especially the miR- 1 5a/miR- 1 6/miRI 95/miR-497 family, miR- 143/miR- 145 and the miR- 1 06-363 cluster appear to be 
important within the known HPV pathogenesis. 

conclusion: This study adds new knowledge to the known pathogenic pathways of HPV and substantiates the oncogenic role of 
HPV in subsets of HNSCCs. 
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Head and neck cancers are the sixth most common cancer in the 
world (Warnakulasuriya, 2009). They consist of a heterogeneous 
group of neoplasms, derived from the epithelium in the upper 
aerodigestive tract. The most common type is squamous cell 
carcinoma (SCC) occurring in the oral cavity, pharynx or larynx. 
Tobacco and alcohol are the predominant risk factors in head 
and neck SCC (HNSCC), but human papilloma virus (HPV) is 
also established to be associated with HNSCC, especially in 
the case of oropharyngeal SCC, that is, in the tonsils and base 
of the tongue. Although the incidence of tobacco- and alcohol- 
induced cancers is declining, there are solid, indications that the 
incidence as well as the prevalence of HPV-induced HNSCCs is 
increasing (Nasman et al y 2009). In Denmark, the incidence of 
oropharyngeal cancers has doubled from 1999 to 2008, in both 
men and women (the Danish Cancer Registry http://www.sst.dk). 
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Human papilloma virus + HNSCC is distinct from HPV-negative 
(HPV — ) HNSCC in a number of ways. Patients with HPV + 
HNSCC tend to be younger and have a lower intake of tobacco 
and alcohol. Moreover, HPV+ tumours are usually smaller 
in T-size but with large cystic nodal involvement (N + ). The 
tumours are frequently non-keratinised (Lajer and Buchwald, 
2010), but normally have favourable prognosis compared with 
HPV- HNSCC (Hafkamp et al, 2008; Lassen et al, 2009). Today, 
it is well established that particular types of HPV cause cervical 
and anogenital cancer. More than 99% of cervical cancer biopsies 
contain high-risk HPV genomes A (Walboomers et al y 1999), 
and cervical cancer is preceded by precancerous lesions, which 
are effectively prevented by vaccination against HPV (Majewski 
et aU 2009). The correlation between HNSCC and HPV is less 
clear. Even though HPV+ cancers are recognised as a distinct 
subgroup within HNSCC (Syrjanen et al, 1983a, b), precancerous 
lesions are not recognised, in oropharyngeal cancers, and the 
role of HPV in HNSCC is still not fully established. The issue 
has major implications for the clinical handling of HNSCC, as 
a causal function of HPV would implicate, that HPV vaccines 
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could be effective in preventing HNSCC in the same way as 
cervical cancers. 

One way to expand the understanding of HPV in HNSCC is 
to define the molecular perturbations in HPV+ tumours and 
compare them to those occurring in cervical cancers. 

The course of malignant transformation in cells infected with 
oncogenic HPV types is incompletely understood, but some 
components are known (Wise-Draper and Wells, 2008). Virally 
encoded E6 binds simultaneously to the ubiquitin/protein ligase 
E6AP and p53, resulting in ubiquitination of p53 and its 
subsequent proteolytic degradation. Human papilloma virus E7 
protein moreover binds and destabilises retinoblastoma protein 
(pRb) and reduces binding to E2F and thereby promotes cell cycle 
progression. Inactivation of pRb also results in a reciprocal 
overexpression of pl6 tumour suppressor protein pl6INK4A, 
which inhibits phosphorylation of the Rb-E2F complex (Ghittoni 
et al, 2010;Pim and Banks, 2010). 

Among the large number of different biomarkers microRNAs 
(miRNAs) have previously been reported to be intimately 
connected to cellular transformation. microRNAs are small 
21-22n non-protein-coding RNAs that regulate mRNA translation 
and decay. microRNAs are involved in terminal differentiation 
processes and have a significant role in the development of human 
cancer. microRNA patterns are tissue-specific and it is normally 
possible to distinguish carcinomas from normal cells. At the same 
time miRNAs are potential therapeutic targets or therapeutic 
agents in future cancer treatments (Pan et al, 2010). We only have 
limited information about miRNAs in HPV+ HNSCC cancer. 
Wald et al (2011) recently examined cultured HPV+ and HPV- 
SCC cells and identified a number of perturbed miRNAs in the 
HPV+ cells. The expression of the HPV- 16 E6 oncogene was 
sufficient to alter miRNA expression (Wald et al, 2011), indicating 
that HPV was directly connected to the observed changes. In 
agreement with these observations, we previously reported a 
difference in the miRNA profiles between the HPV+ and HPV — 
pharyngeal SCCs (PSCCs; Lajer et al, 2011). 

To further characterise these differences, we expand our 
previous observations by applying laser microdissection on 
tonsillar SCC (TSCC) to isolate cancerous and non-cancerous 
tissue and compare the miRNA signatures from HPV+ and 
HPV— tumours with those of CSCCs. The aim of this study is to 
substantiate the distinct oncogenic role of HPV in a subset of 
HNSCC by showing that the miRNA expression in HPV+ and 
HPV— HNSCC is significantly different and that miRNA profiles 
inHPV+ HNSCC show similarities with HPV + CSCC. By making 
a link to CSCC, we aim to explore the common virus-induced 
molecular changes in the cell that leads to the malignant state. 



MATERIALS AND METHODS 
Patients and tumour samples 

All patients with pharyngeal cancer and controls (fresh frozen 
samples as well as formalin-fixed paraffin-embedded (FFPE) 
samples) in this study had surgery performed at the Department 
of Oto-rhino-laryngology, Head & Neck Surgery and Audiology, 
Rigshospitalet, University of Copenhagen. The department is 
responsible for the surgical treatment of head and neck cancer 
patients in East Denmark, the Faeroe Islands and Greenland and 
receives patients from all of Denmark for treatment at highly 
specialised level. Cervical SCC patients and controls had primary 
surgery done at the Department of Gynaecology, Rigshospitalet. 
All fresh frozen samples were taken during primary surgery from 
vital tumour tissue. Patients were included at the consultation 
primary to surgery. 

The study was approved by the Regional Scientific Ethical 
Committee (H-C-2008-080). 



Paraffin- embedded tonsillar samples 

Formalin-fixed paraffin-embedded tonsillar cancers from the 
period 2008-2010 were collected and pl6 immunohistochemistry 
was performed on all samples. On the basis of the amount of 
tumour tissue in the sections, 31 tumour samples were selected for 
laser microdissection and 13 samples from patients with non- 
neoplastic tonsillar disease were included as controls. The latter 
were typically derived from tonsillectomies performed on chronic 
tonsillitis patients. All specimens were reviewed by a head and 
neck pathologist to ensure the pathological diagnosis. 

Fresh frozen oral cavity and pharyngeal cancer 

Fresh frozen tumour samples from the oral cavity and pharynx is 
previously described in Lajer et al (2011). In the previous study, 
PSCC samples had some degree of heterogeneity, probably because 
of the presence of normal tissue in the PSCC biopsies. In this 
study, we therefore included only samples from large tumours that 
were primarily excised and thereby excluded four samples that 
were suspected to contain a high degree of normal tissue. 

Fresh frozen cervical samples 

All cervical cancer patients provided informed consent for the use 
of tissue harvested at surgery. 

A total of 16 patients with cervical cancer and 16 patients with 
non-neoplastic diseases in the uterus had biopsies taken during 
surgery. The biopsies were placed directly in RNA-later and frozen 
within an hour. Samples from 10 HPV+ CSCC patients and 10 
HPV— control patients were used for miRNA profiling. Five of the 
originally 16 cervical cancer samples were adenocarcinomas and 1 
patient sample was very small and these were excluded for the 
miRNA profiling. One control sample was tested positive for HPV. 
This sample and additionally 5 other control samples were likewise 
excluded to obtain the same number of cancer and control samples. 
Information of diagnosis, smoking and alcohol abuse, differentiation 
grade and TNM stage were collected from the patient's charts. 

Microdissection 

The laser microdissection was performed on a PALM MicroBeam 
(Carl Zeiss Microimaging GmbH, Jena, Germany). Each tissue 
block was cut according to the standard protocol into 10 ^m 
sections and placed on MMI RNAse-free frame PEN membrane 
slides. Sections were deparaffinised and stained with cresyl violet 
using the Ambions LCM staining kit made for RNA isolation 
(Ambion, Life Technologies Corporation, Carlsbad, CA, USA). Our 
cutoff level in cancer sections was achieving ^ 80% tumour cells in 
the dissected areas. In the control samples, tissue areas of stratified 
squamous cell epithelium were excised. Tubes were immediately 
frozen at — 80 °C until further processing. During microdissection 
one to four pictures were taken of each cancer and controls for 
documentation (Supplementary Figure 1). 

RNA extraction 

From fresh frozen biopsies, total RNA was extracted and isolated by 
homogenising the frozen samples in TRIzol (Invitrogen, Life 
Technologies Corporation, Carlsbad, CA, USA) using a homogeniser 
(Tissuelyser, Qiagen GmbH, BY, Germany). Briefly, the total RNA was 
extracted from tissue lysate by phase separation with l-Bromo-3- 
chloropropane and RNA precipitation with isopropanol. After 
washing with alcohol the RNA was eluted in RNAse-free water. Total 
RNA concentration of samples was measured using a Nanodrop ND- 
1000 Spectrophotometer (Wilmington, DE, USA). 

In the case of FFPE specimens RNA was isolated using a 
combination of Ambions Recover All and Ambions RNAqueous. 
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DNA extraction 

DNA from fresh frozen biopsies was extracted from the interphase 
of the Trizol and precipitated with ethanol. It was dissolved in 
lysate buffer and proteinase K, and after a second precipitation 
with ethanol, the solutions were transferred to spin columns. After 
washing, DNA was eluted with water. DNA from FFPE specimens was 
isolated from four 10 /mi sections using the Qiagens QIAamp DNA 
FFPE Tissue Kit to spin columns (NucleoSpin Tissue). 

HPV testing 

Immun ©histochemistry for pl6 INK4A was performed on FFPE 
tissue slides pre-treated in DAKO (Dako Denmark A/S, Glostrup, 
Denmark) PT Link, incubated with pl6 (JC8) antibody, Santa Cruz, 
(Santa Cruz Biotechnology Inc., Santa Cruz, CA, USA) and 
visualised using the DAKO Envision Flex + system. Slides were 
categorised as positive, if there was a strong diffuse nuclear and 
cytoplasmic staining of >50% of the tumour cells. PCR for HPV 
was performed with general primers MY09 and MY 11 and GP5 + / 
GP6 + . Validation of the DNA quality and the efficacy of the PCR 
were done with PCR for the house keeping gene GAPDH. Tissue 
with a positive GAPDH PCR was considered HPV+ by PCR if 
either one or two positive primer sets were positive. Cancers were 
classified as HPV+ based on both positive PCR and positive pl6 
staining. Cancers with a negative PCR were considered HPV — 
without regards of the pl6 immunohistochemistry. None of the 
PCR+ samples were pl6 negative. According to these preliminary 
investigations, 10 HPV+ and 10 HPV— cancers were selected for 
laser microdissection. The 10 control tonsils were selected for 
microdissection based on a positive GAPDH PCR. Nine controls 
were HPV— and one was HPV+. All HPV+ samples were 
additionally tested positive with HPV-type 16- and -type 11- 
specific primers. All HPV+ PSCC and FFPE TSCC were positive 
for the type 16-specific primers and additionally negative for type 
11 -specific primers. The 10 HPV+ CSCC samples were typed for 
HPV-type 16 and 18. Seven samples were HPV-type 16 positive 
and three were type 18 positive. 

MiRNA array and data deduction 

A 500 ng of total RNA was labelled with FlashTag Biotin RNA 
Labelling Kits from Genisphere (Genisphere LLC, Hatfield, PA, USA) 
and hybridised to Affymetrix GeneChip miRNA Array according to 
the manufacturer's instructions (Affymetrix, Santa Clara, CA, USA). 

The Affymetrix miRNA Array comprise probes for 847 human 
miRNAs. Four copies of each miRNA probe are distributed on the 
array. To minimise batch variation, each batch was run with equal 
numbers of HPV+, HPV— cancers and controls. Raw data files 
were imported into Affymetrixs miRNA QCTool and normalised 
using the quantiles normalisation and median Polish summarisa- 
tion following a background correction, which corrects for the GC 
content of each particular probe. Log2 intensities of the 847 human 
miRNAs were imported into the Data Analysis software package 
Qlucore Omnics Explorer 2.1 (Qlucore AB 2008-2011, Ideon 
Science Park, Lund, Sweden). Class comparison analysis was 
performed using the Student's f-test. Comparisons of cancers vs 
controls were performed for each of the categories. The principal 
component analysis (PCA) visualisation of the clustering of 
samples using the genes selected in the class comparison was 
performed using the build-in PCA tool in Qlucore Omnics 
Explorer 2.1. 

The number of miRNAs found to be differentially expressed in 
both the HPV+ HNSCC tissues and in the cervical tissue was 
compared with the number of miRNAs differentially expressed in 
both the HPV— HNSCC tissues and in the cervical tissue. To assess 
if the difference was significant, a permutation test was developed 
using the software R (R Development Core Team, 2010). 



qPCR 

qPCR validation was carried out on miRNAs selected from the 
miRNA microarray analysis. cDNA was prepared from 50 to 100 ng 
total RNA from fresh frozen samples and FFPE samples using the 
TaqMan MicroRNA Reverse Transcription Kit (Applied Biosys- 
tems, Life Technologies Corporation, Carlsbad, CA, USA) and 
applied to TaqMan MicroRNA Assays containing predesigned 
primers for mir-16b, mir-363, mir-145,mir-497, mir-381, mir-574- 
3p, mir-20b and mir-15a. Hsa-mir-191 and RNU-44 were used for 
endogenous control. Quantitative RT-PCR reaction was performed 
using the TaqMan Universal PCR Master Mix No AmpEras UNG, 
according to the manufactures instructions, all from Applied 
Biosystems (Life Technologies Corporation). Each amplification 
reaction was performed in triplicate, and median value of the 
three-cycle threshold was used for further analysis. For calcula- 
tions of fold changes we used the 2 ~ Act method. 



RESULTS 

Patient characteristics 

Patient characteristics for the tonsillar patients are summarised in 
Table 1. It was not possible to obtain information about one 
HPV— TSCC patient because the patient was from Greenland and 
did not have a Danish chart. The mean age of all TSCC patients was 
55. The HPV+ cancer patient group and the control group were 8 

Table I Patient characteristics in laser microdissected FFPE tonsillar 
samples 



HPV positive 
TSCC 



HPV negative 
TSCC 



Controls 



Sex 
Male 
Female 

Age (Mean) 

Stage 



IVA 
IVB 
IVC 

Unknown 

Differentiation 
Poorly 
Intermediate 
Highly 
Unknown 

Alcohol abuse 
Yes 
No 

Unknown 

Tobacco abuse 
Current 
Former 
Never 

P/6 
Positive 
Negative 

PCR 
Positive 
Negative 



9 
I 

56 



10 
0 



10 
0 



10 
63 



46 



0 
10 



Abbreviations: FFPE = formalin fixed paraffin embedded; HPV = human papilloma 
virus; TSCC = tonsillar squamous cell carcinoma. 
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Table 2 A total of 36 differentially expressed miRNAs comparing 
HPV+TSCC and HPV- TSCC (P<0.0 1 , q<0. 1 5 and fold changes > 1.5) 



HPV+/HPV- 


P-value 


q- Value 


Fold change 


hsa-miR-363_st 


I.82E-05 


0.007712241 


10 


hsa-miR-2 l_st 


0.001653245 


0.063649926 


5.3 


hsa-miR-l50_st 


0.002013624 


0.068221594 


3.8 


hsa-miR-l46b-5p_st 


0.00329 1 349 


0.092955818 


3.5 


hsa-miR- 1 5a_st 


0.004978216 


0.1 17126349 


3.4 


hsa-miR-20b_st 


0.006161643 


0.13381825 


3.1 


hsa-miR- I46a_st 


0.003292414 


0.092955818 


3 


hsa-let-7g_st 


0.000378428 


0.028654156 


2.9 


hsa-miR-625_st 


0.005353066 


0. 1 2254 1819 


2.8 


hsa-miR- 1 5 5_st 


0.00081 144 


0.036173155 


2.6 


hsa-miR- 1 5 b_st 


0.000624782 


0.03 1 1 28835 


2.6 


hsa-miR-29a_st 


0.000542066 


0.0302208 


2.6 


hsa-let-7f_st 


0.000405476 


0.028654156 


2.6 


hsa-miR- I25b_st 


0.000519276 


0.0302208 


2.2 


hsa-miR-26b_st 


0.003800803 


0.100602505 


2 


hsa-miR-342-3p_st 


0.00475 1 1 8 


0.1 14978563 


2 


hsa-miR- 1 6_st 


0.00010175 1 


0.01231 1862 


1.9 


hsa-miR-768-3p_st 


0.002282882 


0.074369259 


1.8 


hsa-miR-34a-star_st 


0.00643 1 692 


0. 1 36 1 9 1 078 


1.7 


hsa-miR-596_st 


0.000405962 


0.028654156 


1.6 


hsa-miR-20b-star_st 


0.005589569 


0.124588547 


1.5 


hsa-miR-598_st 


6.79E-05 


0.009661928 


1.5 


hsa-miR-200c_st 


0.001818879 


0.066982197 


- 1.6 


hsa-let-7e_st 


0.004208445 


0.104839803 


- 1.9 


hsa-miR- I275_st 


6.0IE-05 


0.009661928 


-2.0 


hsa-miR-l8lb_st 


0.002641257 


0.079898032 


-2.0 


hsa-miR-324-5p_st 


6.84E - 05 


0.009661928 


-2.1 


n sa- m i K-4 z 5 - j p_st 


U.UU 1 7 DADO J 


U.UoozlH- 1 SI 


— 2. 1 


hsa-miR- 1 25a-3p_st 


3.03E-05 


0.008556386 


-2.2 


hsa-miR-99b_st 


0.001 131592 


0.045640877 


-2.4 


hsa-miR-744_st 


0.00021728 


0.023004549 


-2.7 


hsa-miR-877_st 


0.000757677 


0.035652894 


-2.8 


hsa-miR-99b-star_st 


7.45E-06 


0.006308449 


-3.0 


hsa-miR- 1 1 80_st 


0.000541339 


0.0302208 


-3.4 


hsa-miR-3 l_st 


0.003962398 


0.101701547 


-4.7 


hsa-miR- I93b-star_st 


0.00096729 1 


0.040964776 


-5.1 



Abbreviations: HPV = human papillomavirus; miRNA = microRNA; TSCC = tonsillar 
squamous cell carcinoma. 

and 17 years (mean) younger than the HPV— cancer group, 
respectively. Moreover, there was a higher proportion of females in 
the control group compared with the cancer groups. Only one 
female was present in the cancer groups presenting with a HPV + 
cancer. 

The mean age for the CSCC patients was 48 years whereas for the 
cervical controls it was 43. Cancer stages were distributed with two 
patients having stage lb, one patient stage Ha, two patients stage 
lib and five patients stage Illb. Tumour differentiation was three 
patients having intermediate differentiation, four patients having 
low differentiation and three patients were unreported. 

microRNA expression in FFPE tonsillar specimens 

To obtain an overview of the differences and similarities of the 
miRNA expression among the HPV+/— TSCC and normal 
tissues groups, we first compared the miRNA expression of all 
samples in a PCA including all variables (Figure 1). In the 
uncensored PCA plot, samples grouped clearly into three distinct 
groups, demonstrating that miRNA signatures can distinguish 
normal tonsillar mucosa from the malignant tumours. Furthermore, 
miRNA signatures can also separate HPV+ tumours from HPV — 
tumours. In a multigroup comparison analysis with a significance 
level of P<0.01 (false discovery rate (g)<0.06) - 144 miRNAs were 
differentially expressed between the three groups (Figure 1). 

As the HPV+ and HPV— tumours were clearly different, 
we subsequently examined HPV+ and HPV— tumour samples 



separately and compared them with the normal samples in a 
censored class comparison analysis. Using a significance level 
of P-values<0.01 (g<0.05), we found 116 miRNAs to be 
differentially expressed between HPV+ TSCC and controls. 
A total of 111 miRNAs were differentially expressed between 
HPV- TSCC and controls. In a comparison of HPV+ TSCC 
to HPV— TSCC, we found 36 differentially expressed miRNAs 
(P-value < 0.01, g-value<0.15 and fold changes > 1.5). In all, 22 
miRNAs were upregulated and 14 were downregulated (Table 2). 

microRNA expression in cervical specimens 

Differentially expressed miRNAs from the cervical specimens were 
visualised in a PCA plot, including all variables, and it was noticed 
that the CSCC samples and control samples clustered into two 
distinct groups, demonstrating that the miRNA expression 
differentiates normal cervical tissue and CSCC (Figure 2). In a 
censored two group class comparison analysis, using a significance 
level of P-values<0.01 (g<0.05), 148 miRNAs were differentially 
expressed between the CSCC and cervical controls (Figure 2). As 
CSCC is known to be driven by HPV, we used the 148 miRNAs 
differentially expressed in CSCC as a reference for HPV-mediated 
changes (Supplementary Table 1). 

Eight miRNAs were selected for validation by qPCR. MiR-363, 
miR-16, mir-574-3p, mir-20b and mir-15a were used for the FFPE 
samples. Whereas the fresh frozen samples (HPV+ PSCC, HPV — 
oral SCC (OSCC), pharyngeal/oral controls and all cervical SCC and 
controls) were additionally tested for mir-145, miR-497 and mir- 
381. In each case, cancer vs controls was compared. The results of 
the qPCR validation are shown in Supplementary Figure 2. By qPCR 
miR-363 was also found upregulated in the CSCC samples and miR- 
497 was also found downregulated by qPCR in the HPV- PSCC. In 
all other cases qPCR results confirmed the array data. 

Comparison of all HPV + and HPV- HNSCCs with CSCC 

In an attempt to define common HPV-controlled miRNAs, that are 
independent from tissue type and anatomical localisation, we 
compared the differentially expressed miRNAs of HPV+/ HPV — 
HNSCC and their respective controls and compared them to the 
miRNAs in the HPV+ CSCCs. The list of differentially expressed 
miRNAs from the HPV+ FFPE TSCC was merged with our 
previously published differentially expressed miRNAs from fresh 
frozen HPV+ PSCC and defined as HPV+ HNSCC miRNAs. In 
this group, we identified 144 differentially expressed miRNAs. 
Likewise, the list of differentially expressed miRNAs of HPV — 
FFPE TSCC was merged with previously published differentially 
expressed miRNAs in fresh frozen HPV— PSCC and defined as 
HPV- HNSCC miRNAs. Here, we identified 140 differentially 
expressed miRNAs. Subsequently we compared the miRNA 
expression profiles in the HPV+ and HPV— HNSCC group with 
the miRNA expression profiles in CSCC. 

How these two sets of miRNAs overlap with each other and the 
miRNAs identified in the cervical samples is illustrated in the Venn 
diagram (Figure 3). 

Of particular interest is the number of miRNAs shared 
exclusively by HPV+ HNSCC and CSCC compared with the 
number of miRNAs shared exclusively by HPV— HNSCC and 
CSCC. As seen from the diagram d=n + —n_ =23-6 = 17. The 
distribution of d under the null hypothesis that HPV+ and 
HPV— samples are exchangeable is shown in Supplementary 
Figure 3. In 4 out of 1000 permutations, the value of d turned out 
> 17, giving the test a P-value of 0.004. The permutation test was 
also run with a harder cutoff (false discovery rate q < 0.05 by the 
Benjamin-Hochberg adjustment for multiple testing) with a 
similar result (P< 0.032). 

To be absolutely certain of identifying miRNAs that are 
exclusively present in HPV+ tumours, additionally 7 miRNAs 
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Figure I (A) Principal component analysis (PCA) plot including all variables. Human papilloma virus— TSCC (yellow), HPV+ TSCC (purple) and 
tonsillar controls (black). (B) Principal component analysis plot, visualising differentially expressed miRNAs in HPV+ , HPV— and tonsillar control, subset of 
variables. (P<0.0l, q<0.06). (C) Heat map: three distinct groups were observed. Human papilloma virus — TSCC (yellow), HPV+ TSCC (purple) and 
tonsillar controls (black). The heatmap and PCA visualisation was performed using the Qlucore Gene Expression Explorer (http://www.qlucore.com). 




were excluded from the 23 miRNAs identified only in HPV+ 
HNSCC and CSCC, as they have previously been identified in 
HPV— OSCCs by our study group (Lajer et al, 2011). microRNAs 
denoted with a star was likewise excluded. That leaves 14 miRNAs 
to be defined as HPV core miRNAs, which are listed in Table 3. In 
all, 4 were upregulated and 10 were downregulated. The 
upregulated HPV core miRNAs were miR-lOb, miR-15a, mir-20b 
and miR-16 and among the downregulated HPV core miRNAs 
there were miR-139-3p, miR-145, miR-381 and miR-574-3p. 

microRNAs that were common for all SCC independent of HPV 
status and compared with their respective normal tissue could 
therefore be defined as SCC cancer-specific miRNA. These are 
listed in Table 4 and include the miRNAs that were identified in 
HPV+ HNSCC/HPV- HNSCC and CSCC as well as previously 
identified in OSCC (Lajer et al y 2011). Finally, 45 (from the Venn 
Diagram) miRNAs were common between the HPV+ HNSCC and 
HPV- HNSCC but not differentially expressed in the CSCC. These 



miRNAs seem to be cancer-specific miRNAs that are dependent on 
an anatomical localisation within the head and neck. 



DISCUSSION 

The presence of HPV defines a distinct group of head and neck 
cancers exhibiting epidemiological, histopathological, clinical and 
prognostic differences compared with the typical HNSCC. 
Although HPV is acknowledged as the cause of the vast majority 
of cervical cancers, it is still debated which role HPV has in 
HNSCC cancers. To exploit the role of HPV in HNSCC, it is 
important to characterise the molecular basis of HPV+ and 
HPV — HNSCC which can be further strengthened by comparison 
with the carcinogenic mechanisms to HPV+ CSCC. 

This study is the first to make a direct correlation of miRNA 
expression in HPV+ HNSCC to CSCC in the same settings and in 
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Figure 2 (A) Principal component analysis plot including all variables. (B) Principal component analysis plot, visualising differentially expressed miRNAs in 
CSCC compared with cervical controls. (C) Heat map (P<0.0l, q<0.05). Cervical squamous cell carcinoma (CSCC) (purple), cervical controls (black). 



that way, it demonstrates the impact of HPV on the miRNA expres- 
sion in solid malignant tumours irrespective of the anatomical site. 

The demonstration of HPV-dependent perturbations on miRNA 
expression adds to the complexity of HPV-induced transformation. 
Human papilloma virus elicits widespread changes in cellular gene 
expression, which may at least partially be responsible for the 
distinct clinical behaviour of HPV+ cancers. 

In this study, we show that HPV+HNSCC have a distinct 
miRNA profile that clearly distinguishes them from HPV — 
HNSCC. Furthermore, there are several similarities in the miRNA 
expressions of SCC in spite of two different anatomical localisa- 
tions. microRNAs are highly tissue specific (Rosenfeld et al, 2008; 
Rosenwald et al, 2010). As expected we found a high degree of 
concordance between HPV+ and HPV— HNSCC and a much 
lesser concordance between CSCC and HNSCC. Interestingly, 
however, the permutation test confirmed that HPV+ HNSCC and 
HPV+ CSCC were significantly (P = 0.004) more alike compared 
with HPV- HNSCC and CSSC. 



It was striking that the miRNAs identified as HPV core miRNAs 
appear to involve closely related miRNAs united in a few clusters 
and miR- families and they have distinct interactions within theE6 
p53 and E7/pRb pathways of HPV. Figure 4 illustrates the 
proposed interactions of the HPV core miRNAs in the HPV- 
induced malignant pathogenesis. 

Among the HPV core miRNAs, we identified five members of 
the miR-family c miR-15a, miR-16, miR195, miR-424 and miR-497\ 
We observed increased levels of miR- 15a and miR-16 and 
decreased levels of miR- 195 and miR-497 (the latter two miRNAs 
were also significantly reduced in HPV— HNSCC although to a 
lesser degree). These miRNAs target mRNAs-encoding proteins 
involved in Gl/S checkpoint and they have been considered as 
tumour suppressors. Studies show that cell lines lacking functional 
Rb (as in HPV-driven cancer) are resistant to mir-15a-16 cell cycle 
arrest and reintroduction of functional Rb resensitises the cells 
(Bandi et al, 2009). Mir- 15a and mir-16 are often downregulated in 
cancer and their expression is inversely correlated with the 
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Figure 3 Venn diagram illustrating the number of miRNAs differentially 
expressed in the three groups (HPV+ HNSCC, HPV- HNSCC and 
CSCC) compared with their normal counterpart and how they overlap 
with each other (P<0.0l). 



Table 3 Human papilloma virus core miRNAs 
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Abbreviations: CSCC= cervical squamous cell carcinoma; HNSCC = head and neck 
SCC; HPV = human papilloma virus; miRNA = microRNA; OSCC = oral squamous 
cell carcinoma; TSCC = tonsillar SCC. P<0.0 1 . MiRNAs with star or formerly shown 
to be present in OSCC is removed. 



expression of cyclin Dl. In non-small cell lung cancer cell lines the 
cyclines Dl, D2 and El have been shown to be directly regulated by 
the level of mir-15a/mirl6 and overexpression of these induces cell 
cycle arrest in the Gx-Gq phase. 

The most markedly and consistently upregulated miRNA in 
HPV+ TSCC was miR-363 (Supplementary Table 3). This is in 
concordance with both our previous study and the study by (Wald 
et al, 2011), demonstrating a considerable HPV 16 E6-dependent 
upregulation of miR-363 in HPV+ cell lines. In agreement with 
other studies (Wang et al, 2008; Pereira et al, 2010), miR-363 was 
not found differentially expressed in our CSCC specimens in the 
array data, though it was upregulated by qPCR. The miR- 106-363 
cluster was recently reported implicated in T-cell leukaemia 
(Landais et al, 2007), and overexpression of the miR- 106-363 
cluster induced anchorage independent growth of NIH/3T3 cells. 
Five of the six miRNAs in this cluster were upregulated in the 
CSCC samples (miR- 106a, miR 18b, miR-20b miR92a and miR-92b; 



Table 4 Cancer-specific miRNAs, that is, differentially expressed 
miRNAs in cancer vs controls in all three categories (HPV+ HNSCC 
(FFPE TSCC/FF PSCC), HPV- HNSCC (FF PSCC/FFPE TSCC) and CSCC 
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Abbreviations: CSCC = 


cervical squamous cell 
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paraffin embedded; HNSCC = head and neck SCC; HPV = human papilloma virus; 
miRNA = microRNA; OSCC = oral squamous cell carcinoma; PSCC = pharyngeal 
SCC; TSCC = tonsillar SCC. formerly identified in HPV- OSCC. 

Supplementary Table 1). Besides miR-363, also miR-20b was 
upregulated in the HPV+ HNSCC. 

We identified miR- 145 in the core set of HPV-controlled 
miRNAs as it was downregulated in CSCC as well as in HPV + 
HNSCC. MiR-143 was also downregulated in both HPV+ HNSCC 
and in CSCC, but did not reach a significance level of P<0.01 
(P = 0.011) in the cervical samples. MiR-143 and miR-145 have 
previously been shown to suppress cell growth in CSCC and they 
have consequently been implicated in cervical carcinogenesis 
(Wang et al, 2008). Repression of miR- 143- 145 could be the result 
of the HPV oncogene E6 as shown in Figure 4. Even though a 
reduction of miR- 143 and miR-145 seems to be related to HPV in 
head and neck cancers, reduced expression of these miRNAs have 
been described in other tumour types unrelated to HPV presence 
(Michael et al, 2003; Akao et al, 2007a, b). 

Reduced levels of miR- 195, miR- 143 and miR-145 may be 
secondary to E6-dependent loss of p53/p63/p73 that normally 
regulate these miRNAs (Boominathan, 2010). Zheng and Wang 
(2011)also describes how the viral oncoproteins E6 and E7 
deregulate the expression of the miR- 15/ 16 cluster, the miR- 17-92 
family and miR-106b/93/25 cluster via the E6-p53 and E7-pRb 
pathways in HPV cervical cancers (Zheng and Wang, 2011). This 
study supports these findings in HPV+ HNSCC as well. 

Although HPV+ and HPV— tumours could be clearly 
distinguished, the differences in miRNA expression were diminu- 
tive. In general, the presence of HPV DNA in the cancer samples 
was associated with downregulation of miRNAs in CSCC and 
HPV+ HNSCC (10 of 14, Table 4). 

A difference in gene expression profiles between HPV+ and 
HPV- HNSCC as well as a resemblance of the HPV+ HNSCC with 
CSCC have been reported in another study (Pyeon et al, 2007). We 
can now deduce that both miRNA and mRNA expression profiles 
clearly display differences in HPV+ and HPV- HNSCC. In 
addition, the miRNA and mRNA expression in HP V+ HNSCC 
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Figure 4 The schematic representation of the relation between HPV-regulated miRNAs and factors involved in malignant transformation caused by the 
viral oncoproteins E6 and E7. Putative miRNA-regulated proteins and miRNA controlled by transforming factors are indicated in grey and miRNAs are 
shown in blue boxes. Binding of E7 to pRb leads to degradation of pRb and release of transcription factor E2F, which drives the cell into S-phase and cell 
division. Interaction between E7 and p2l results in inhibition of the cyclins, which are targeted by miR-15, miR-16, miR-195, miR-497, and this leads to 
inhibition of cell cycle arrest and contribute to carcinogenesis. Ubiquitination by E6 and ubiquitin ligase E6AP (potentially targeted by miR- 1 39-5p) leads to 
p53 degradation and subsequently impaired apoptosis following DNA damage. p53 is also potentially targeted by miR-574-5p and miR-38 1 , and p53 in turns 
regulates transcription of miR-15, miR-16, miR- 1 43, miR- 1 45 and miR-195 By association with E6/E6AP, NFXI is degraded and results in activation of 
hTERT, which leades to cellular immortalisation. E6 degrades the tumour supressor PTPN I 3 and loss of PTPN I 3 further enhances the oncogenic Ras 
signalling. Ras further regulates miR- 1 43 and miR- 1 45 transcription by binding to the promoter. 



exhibit similarities to the known HPV-induced CSCC. This further 
elucigcqias the oncogenic role of HPV in different anatomical sites. 

This study reflects the global change in miRNA expression of the 
different cancer types. The intention is, as suggested by Gallagher 
et al (2010), to investigate the collective changes in many miRNAs 
(global profiling), which may be the most interesting biologically 
parameter to consider - compared with choosing one miRNA for 
further functional analysis in a cancer cell line, which may be far 
from the conditions in vivo. It is most likely that miRNAs work 
cooperatively in vivo and that physiological regulation of a single 
protein by a single miRNA may be a rather rare occurrence. In 
accordance with this, we observed that changes of miRNA 
expression often implicated a cluster of miRNAs or a whole 
miRNA family. 

Most of the HPV core miRNAs identified in this study have not 
been investigated in relation to HNSCC, as most miRNA studies 
regarding HNSCC have not taken the HPV status into considera- 
tion. Further functional studies of these miRNAs in HPV+ and 
HPV- HNSCC cell lines will be needed in the future. 

It is the combination of the HPV core miRNAs that seems to be 
closely related to HPV infection, regardless of the fact that the 
individual miRNAs may not be unique to HPV-driven cancer. 
Consequently, some of the HPV core miRNAs have also been 
identified in other types of HPV— cancers. This aspect will need 
further investigations. 



CONCLUSION 

In this study, we show that HP V+ HNSCC has a distinct miRNA 
profile compared with HPV — HNSCC. Furthermore, miRNA profiles 
in HPV+ HNSCC are significantly more similar to miRNA profiles 
of CSCC than HPV- HNSCC. A set of HPV core miRNAs were 
identified. Of these especially the miR-15a/miR-16/, miR-143/miR-145 
and the miR- 106-363 cluster appear to have targets within the known 
HPV pathogenesis. This adds new knowledge to known pathogenic 
pathways of HPV and substantiates the oncogenic role of HPV 
in subsets of HNSCCs. The identification of miRNAs associated 
with different HPV-mediated cancers is important, as they may be 
useful as general and HPV-specific tumour markers. In addition, 
they may prove useful in future individualised treatment strategies 
and development of new targeted therapies. 
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